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I— I ■ Abstract 

Qh| We revisit the formulation of the principle of minimal flavor violation (MFV) in the 

pi i' minimal supersymmetric extension of the standard model, both at moderate and large 

^ . tan/3, and with or without new CP-violating phases. We introduce a counting rule 

which keeps track of the highly hierarchical structure of the Yukawa matrices. In this 
^ manner, we are able to control systematically which terms can be discarded in the soft 

^ . SUSY breaking part of the Lagrangian. We argue that for the implementation of this 

counting rule, it is convenient to introduce a new basis of matrices in which both the 
squark (and slepton) mass terms as well as the trilinear couplings can be expanded. 



f^ ■ We derive the RGE for the MFV parameters and show that the beta functions also 

r — I , respect the counting rule. For moderate tan/3, we provide explicit analytic solutions 

^^ I of these RGE and illustrate their behaviour by analyzing the neighbourhood (also 

switching on new phases) of the SPS-la benchmark point. We then show that even in 
the case of large tan (3, the RGE remain valid and that the analytic solutions obtained 
for moderate tan/3 still allow us to understand the most important features of the 
^ ' running of the parameters, as illustrated with the help of the SPS-4 benchmark point. 
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1 Introduction 

In the Standard Model, renormalizability restricts the possible sources of flavour violations 
to only one matrix, the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Its almost diagonal 
structure (like the values of all the other free parameters of the model) remains unexplained. 

In extensions of the Standard Model, as soon as new degrees of freedom appear, the 
possibilities to generate transitions among flavours increase very rapidly, and may give rise 
to a richer phenomenology than what the Standard Model alone would allow. The wealth 
of recent experimental results produced at kaon and B factories shows, however, that if new 
degrees of freedom exist just above the electroweak symmetry breaking scale, their influence 
on low-energy flavour physics is smaller than one would naively expect, below the current 
experimental sensitivity. This imposes a nontrivial constraint for model building, and forces 
one to impose some sort of protection against flavour violations. 

A convenient way to do this, without excluding completely the possibility to have new 
effects in flavour physics, is the principle of minimal flavour violation (MFV) [U [2] (see also 
Ref. [3]). According to this, even in extensions of the Standard Model, the only source of 
flavour violations is in the Yukawa matrices, and since one of them can always be diago- 
nalized, in the CKM matrix. In its latest, more complete implementation, the principle is 
formulated as a symmetry: one starts from the observation that the large global flavour 
symmetry group of the Standard Model in the absence of the Yukawa couplings is saved 
even in their presence if they are promoted to the status of spurious, i.e. if they are assigned 
transformation properties under the flavour group. The principle of MFV requires that the 
same symmetry holds even in extensions of the Standard Model - the only allowed spurious 
being the Yukawa matrices. 

In this paper, we concentrate on supersymmetric extensions of the Standard Model with 
minimal fleld content and exact R-parity (MSSM) (for a nice introduction, cf. Ref. [1]) and 
discuss the principle of MFV within this framework. Our main points are: 

1. We observe that, if one considers the Yukawa matrices as dimension-zero spurious, and 
allows any power of them to appear in local operators, imposing MFV on the MSSM 
does not restrict the number of free parameters of the model, but amounts to a mere 
reparametrization. Still, if one requires that the coupling constant of the model are 
of order one, some of them are irrelevant for the phenomenology and can therefore be 
dropped. 

2. In order to decide systematically which terms are irrelevant and which ones should be 
kept, we flnd it convenient to introduce a counting rule, and use as expansion parameter 
A, the Cabibbo angle. We use this expansion parameter to take into account not only 
the highly hierarchical CKM matrix, but also the highly hierarchical quark and lepton 
masses. Each of the MFV parameters will therefore be assigned an order in A - once 
one has set the accuracy of its calculation to a certain level, (9(A"), it is immediate to 
see which MFV parameters should be kept. 

3. Since imposing MFV amounts to a mere reparametrization of the generic MSSM, it 
is obvious that the principle is renormalization group (RG) invariant - what is not 
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guaranteed is that the coefficients will keep their order in A during the running. We 
will rewrite the renormalization group equations (RGE) in MFV form, checking that 
the beta functions also respect the same counting rules as the parameters themselves. 
As such, this is only a necessary, but not yet sufficient condition to prove that the 
MFV principle is truly RGE invariant. We will therefore study the solutions of the 
RGE numerically, and show that indeed they respect the counting rules when evolving 
from the high to the low scale - none of the irrelevant parameters at the high scale 
may become so large during the RG evolution that it becomes of phenomenological 
importance at the low scale. We will also show that the converse is not true: generic 
MFV-like boundary conditions at the low scale will not necessarily evolve to a MFV- 
like MSSM at the high scale. Therefore, the assumption that the MFV hypothesis in 
the MSSM is valid at different scales makes it even more restrictive at the low scale. 

4. We discuss in detail the possible new CP-violating phases allowed by the MFV hy- 
pothesis, and analyze their behaviour under the running. We show that they tend to 
vanish at the low scale - how fast they do that depends on the initial conditions at the 
high scale. 

While this paper was being completed, a preprint appeared which also analyzes the 
behaviour of MFV models under running [5j. There, the numerical analysis is performed 
with the help of SOFTSUSY [6], one of the available codes which allow one to run the 
MSSM with a generic flavour structure according to the full RGE to two loops. In this 
paper, the authors start with MFV-compatible initial conditions at the GUT scale, evolve 
all the parameters down to the electroweak scale, and project back the model on the MFV 
parameters, with the help of a ffi. Their analysis is valid for moderate tan/5, and only for 
real MFV coefficients. In this manner, they find out that the MFV parameters have quasi 
fixed points at the low scale. We will con&m their finding, and also provide an analytical 
explanation for this behaviour. Furher, we will perform the analysis also for large tan f3, and 
in the presence of new CP-violating phases in the MFV expansions. 

2 Revisiting minimal flavour violation 

2.1 Definition of minimal fiavour violation 

Gauge interactions in the Standard Model are fiavour blind. If one sets the Yukawa matrices 
to zero, the Standard Model becomes invariant under a large global symmetry group Gp ~ 



GF = Gq®G,® U{1)b ® U{1)l ® U{1)y ® U{1)pq ® f/(l) 



E 



R 1 



where 



G, = SU{?,)Q^®SU{?,)u^®SU{'i)D^, G, = SU{-i)L^®SU{?,)E^ • (2) 

The five f/(l) factors have been decomposed in the three which remain a symmetry even in 
the presence of Yukawa interactions (related to baryon and lepton number and hypercharge) , 



and the remaining two. Following Ref. [2j, we write these as a phase transformation affecting 
Df{ and E^ at the same time, the Peccei-Quinn symmetry of the two-Higgs doublet model 
(denoted here by U{X)pq) and one affecting only Er. 

The Yukawa matrices break G q ® G ^ ® U {1) pq ® U {1) Ej^'- 

Cy = Ur Y„ QlH + Dr Yrf QlH, + Er Ye LlH, + h.c. , (3) 

where He = iT2H*. As observed in Ref. [2|, the Standard Model remains formally invariant 
under Gq ® Ge in the presence of the Yukawa matrices if these are promoted to spurion fields 
transforming as 

Y„~(3,3,l), Y,~(3,l,3) under G,, (4) 

and 

Ye ~ (3, 3) under Gi . (5) 

The symmetry is broken whenever the Yukawa matrices are frozen at a certain value - on 
the other hand, different forms of the Yukawa matrices which are related by Gq and Gi 
transformations are physically equivalent. In what follows, we will choose the following 
background values 

Y„ = XuV , Yd = \d , Yg = Ae , (6) 

where A„ = dmg{yu, Vc, Vt) , ^d = diag{yd,ys,yb) and Ae = diag{ye,y^,yr), and V is the CKM 
matrix. 

An extension of the Standard Model is said to respect MFV if it is symmetric under 
Gq ® Gg in the presence of Yukawa spurious. While new matter fields in such an extension 
are of course allowed, one is not supposed to introduce new spurion fields beyond the Yukawa 
matrices. 

2.2 MFV in the MSSM: a reparametrization of the soft SUSY 
breaking terms 

In a super symmetric extension of the Standard Model, the superpotential automatically 
satisfy the MFV principle. On the other hand, MFV strongly constrains the soft supersym- 
metry breaking terms. We will illustrate this statement by considering the mass term for 
the left-handed squarks 

C^.^ ^ -gtm^ . g , (7) 

and showing first that MFV amounts in principle to a reparametrization of a generic hermi- 
tian 3x3 matrix, and later that if one excludes the possibility of having enormous coupling 
constants, MFV is indeed quite constraining. MFV requires this term to become formally in- 
variant under Gq, hence it must transform like (8, 1, 1). Since we are not allowed to introduce 
new spurious, we have to obtain this transformation property with the help of the Yukawa 
matrices, as with terms like Y^Y^ or Y^Y^. Moreover, one can construct invariants under 
Gq also with the help of e-tensors. Such terms have recently been systematically studied 
in Ref. [H], and permit to extend the MFV principle to the R-parity violating interactions. 
Interestingly, MFV alone is then sufficient to prevent the proton from decaying too rapidly 
- a fact which lends additional support to the validity of the MFV hypothesis at low energy. 



For the R-parity conserving sector of the MSSM, on which we concentrate in the present 
paper, these e-tensor terms have been shown to be very suppressed [H] . Hence, ignoring them 
altogether, MFV permits to write niQ as the infinite sum: 

m^ = ai 1 + 61 Yt Y„ + 62 yJy^ + ci Yt Y„ Yt Y, + C2 YjY.YjYrf + . . . (8) 

Actually, the infinite sum collapses on its first few terms: the Yukawa couplings are 3x3 
matrices, hence they respect the corresponding Cayley-Hamilton identities. The hermitian 
matrix iuq does not contain more than nine independent real parameters and it can be 
shown that the sum in Eq. ([8]) spans the space of hermitian matrices [9] . 
The Cayley-Hamilton identities read for 3x3 matrices 

X^ - (X) X^ + ix ((X)' - (X^)) - det X = , (9) 

and can be rewritten in terms of traces only if the determinant is expressed as 

detX = i(X='>-i(X)(X^> + i(X)^ . (10) 

In other words, all powers of three or more of a combination of Yukawa matrices, X"^^, 
can be eliminated in terms of only X^, X, 1, with coefficients involving the trace of X^ and 
X. Further, identities involving two (or more) different combinations, A and B say, can be 
found by substituting X = aA + 6B in Eq. ([9]) and extracting a given power of a and b. For 
example, a relevant identity is 

A=^B + ABA + BA^ = A^ (B) + (AB + BA) (A) + A ((AB) - (A) (B)) 

+ 1b ((A^) - (A)^) + i (B) ((A)^ - (A^)) + (A^B) - (A) (AB) , (11) 

with A = Yt Y„ and B = Y^Y^. 

Taking into account these identities, the most general expression for iuq respecting MFV 
becomes 

mj = zil + z^YiY^ + z^YlYd + z^iYlY^)^ + z^iY^YdY 
+ z, (Y^Y.Yt Y„ + h.c.) + z^Yt Y.Y^Y.Yt Y, 

+ ZgY^Y^Yt Y„YtY, + z, {{YiY.,f{Y\Y,f + h.c.) 

+ iz^^{Y\YdYiYu - h.c.) + iz^^ ((YtY„)2Y;Yrf - h.c.) 

+ iz,2 ((YjY,)2Yt Y„ - h.c.) + iz,, {{YiY^f{Y\Y,f - h.c.) 

+ iz,, (Yt Y„YtY,(Yt Y„)2 - h.c.) + iz,, (Y^Y^Yt Y.(YtY,)2 - h.c.) 

+ z^i6 (Yt Y.(Y^Y,)2(Yt Y„)2 - h.c.) + iz,, {y\Y .{YIY ^f {Y\Y ,f - h.c.) , (12) 



with the Zi being real parameters. A generic 3x3 hermitian matrix can be described by 
nine real constants - in Eq. (1121) . itiq is expressed in terms of seventeen real constants, so 
that eight of them must be linearly dependent. Even if we do not specify the linear relations 
which allow one to eliminate eight of these constantdj, it is clear that MFV amounts to a 
mere reparametrization of the soft SUSY-breaking terms of the MSSM, since the original 
expansion in Eq. (JHl) contains a basis and the Cayley-Hamilton relations are exact. A similar 
argument can be used also for all other terms. 

What is special about the MFV parametrization is that if all the Zis are of the same 
order of magnitude, the structure of rriQ is highly non-generic. Conversely, if one writes 
down a generic mg matrix and projects it on the MFV basis, the coefficients Zi so obtained 
will typically span many orders of magnitude. We define extensions of the Standard Model 
respecting MFV by the additional requirement that the coefficients appearing in front of the 
various MFV terms are of the same order of magnitude. 

2.3 Counting rules and a new basis 

If one takes all the Zi coefficients to be of the same order of magnitude, several of the terms in 
Eq. flT^ (and in the analogous ones for the other soft SUSY-breaking terms) can be disposed 
of. In this subsection, we will discuss how to do this in a systematic way, and will argue 
that it is more convenient to change basis in order to work with MFV. For example, taking 
into account the actual values of the Yukawa coefficients of the up, charm and top quarks, 
we conclude that the two matrices 

(YtY.)' = VtA^y, YtY„ = VtA2v^ (13) 

are proportional to each other up to a correction of relative order 0{y'^/yf), 

(YtY„)^ - vM-^u = yhXW2, + 0{yt) , (14) 

which is usually neglected since yl <til. A similar argument can be applied to Y^Y^, since 
y1 <til also. We will therefore never include any power of Y^^Y^ or Y^Y^ in our analysis, 
and keep the latter to allow tan/3 to be large (remember that the two Higgs doublets of 
the MSSM separately give mass to the up and down-quarks: VuXu = diag{mu,mc,mt) and 
VdXd = diag{md,ms,mb), with Vu,d the two Higgs vacuum expectation values, and tan/? = 



^Note that such a linear relation can be nontrivial and may involve large coefficients. We stress that 
the Cayley-Hamilton identities do not involve any large numerical coefficients, and so do not upset the 
assumption that the MFV coefficients are of order one. 



m^ 
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Vulvd)- With this approximation, the MFV version of the soft SUSY-breaking terms rea'^^'^' 

ai + 6iYt Y„ + foaY^Yrf + (foaY^Y.Yt Y„ + h.c.) 
a2 + Y„(65 + ciYlY„ + csYJy, + (caYjY.YlY^ + h.c.)) Y^ 
% + ^d[h + C4YtY„ + csYjYrf + (ceYjY.Yl Y„ + h.c.)) yJ 
A^ = AoY„(a4 + &7Y^Y, + cyYtY, + C8Y;fY,YtY„ + C9YtY,Y;fY,) , 
A^ = Ao Y, (as + hYlYu + cioY^Y^ + cn Y;^ Y^Y^ Y„ + ciaY^ Y^Y^Y^) . (15) 

Until now, we have used only the fact that the Yukawa matrices are highly hierarchical 
along their diagonal. They are, however, also hierarchical in their off-diagonal structure, 
and taking this into account leads to further simplifications. In order to do this in a sys- 
tematic way, we use as expansion parameter the Cabibbo angle A = 0.23, which appears in 
the Wolfenstein parametrization of the CKM matrix as follows (to leading order in A - in 
subsequent calculations we will always include higher orders also) 



V 



1 A AX^ (p - i7]) 

-A 1 A\^ 

AA3(1 -p-ir^) -A\^ 1 



(16) 



We then adopt the following counting conventions for the quark mass ratios at the elec- 
troweak scale /i = M^: 



^~0(A^), ^~0(A^), yt-^Oil) 
rrit nit 



;i7) 



rrid 



0{\' 



m.o 



mt mt 

For example, the difference in Eq. 014p . 



Y^Y Y"^Y — ?/^Y^Y 



C(A^ 



mt 



0{X^) . 



C(A8) C(A^) C(A9) 
0{\') C)(A«') C)(A8) 
C(A9) 0(A8) C(AiO) 



represents a correction of at least (9(A^) (in every entry) to each of the two terms, cf. 

YlY„ = I O (X^) O (X^) O (X^) I . (19) 



u ■"- u 



0(A6) C(A5) 0(A3) 
C)(A5) 0{X^) 0(A2) 
C(A3) C(A2) 0(1) 



This shows that while MFV can indeed be viewed as a reparametrization, cf. Eq. 0121) . it is 
on the other hand a very special one, because the basis in the linear space of 3 x 3 hermitian 



^Tlie numbering of the coefRcients follows the choice of Ref. P] whenever possible. The term 64 present 
in that paper has to be equal to 63 in order to satisfy the hermiticity of the matrix hIq. 
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matrices on which ttiq is projected is almost a singular one. Several of the basis vectors 
are almost parallel to each other, and their difference is tiny in comparison to either of the 
vectors. Since our aim here is to link the MFV concept to a counting rule, and define clearly 
which terms should be kept and which should be neglected, it is a lot more convenient 
to use a basis of vectors which are as little U aligned to each other as possible. In this 
way, contributions which are small and may (or may not) be neglected will not have to be 
searched for in small differences between similar contributions, but will be clearly identified 
and separated from the rest. In order to illustrate this concept, we come back to Eq. flT^ 
and observe that the large piece in both Y^Y„ and its square is proportional to the matrix 
V^Vsj, whereas the small one, as indicated in Eq. IHM is proportional to V^V2j. So, instead 
of using Y^jYu and (Y^^Y^)^ as basis vectors, and allow both to have coefficients of order 
one, we find it more convenient to use V^*V3j and V^V2j, and say that the first can have a 
coefficient of order one, but the latter should have it of order y"^ ~ A^, and (depending on 
where one sets its accuracy) can therefore be neglected. 

If we follow the same logic for Y^Y^ and its square, we find that in this case we should 
rather use as basis vectors the matrices 6i36j3 and 6126 j2, and that the coefficients of these 
terms should be of order yl and y"^ respectively (how to translate this into a power of A 
depends on tan/5). Taking into account all the possible structures which can emerge, and 
which can all be multiplied with each other, we are led to consider a set of sixteen matrices, 
which form a closed algebra under multiplication: 



Xi = S^iSsj X5 = 6siV3j Xg = V^i5-ij Xi3 = V^iVsj 

X2 = S2iS2j Xq = 62iV2j Xio = V^*^2j -^14 = V'M^'ij 

X3 = S3iS2j X-j = S3iV2j Xu = V^i52j Xi5 = V^^2j 

X4 = S2i63j Xg = 62iV3j X12 = V^*(53j XiQ = V-^V^j 



(20) 



Notice that all these matrices have at least one entry (almost) equal to one, so that they can 
all be counted as of order one. 

If we now write Eq. (1151) (or even start from the earlier stage given by Eq. (IT2l) ) in this 
basis and assume that all a,, bi and c, coefficients are of order one, we can immediately read 
out the size of the coefficients in front of the Xj matrices and decide which one to keep and 
which not. In particular, if we drop terms of order A^ ~ 10~^ and higher, we can reduce the 
soft SUSY-breaking terms to the following 

niQ = ai + X1X13 + yiXi + 1/2X5 + y^Xg , 

m^ = 0.2 + X2X1 , 

m?) = «3 + VzXx + ti'iXs + w\X4 , 

A^ = S4X5 + 2/4X1 + ti;2X6 , 

A^ = ~ar,Xx + 2/5^5 + W3X2 + W4X4 , (21) 

where, for simplicity, we have absorbed the overall scales vn^ and Aq into the coefficients. 
Since the matrices Xj are all of order one, the coefficients now carry the order in A, and we 



•^Strictly speaking one could - in contrast to what we will do - consider a fully orthogonal basis of matrices 
(like efj = ^\^\)- In this way, however, one would lose track of the fact that the physically relevant degrees 
of freedom of the Yukawa matrices are contained in two diagonal and one unitary matrix. 
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have reflected this in the symbols which identify them. Relative to the leading terms, the 
two Xj's are of order one, a^ (which now incorporates y^ - while 0,4 incorporates yt), as well 
as the ?/j's can become of order one if tan/3 ~ A~'^, and the Wj's are suppressed by at least 
two powers of A, even when tan/5 ~ A~^. More specifically: 

(1\ Ci^ (I4. 

^~0(Ai°4), ^r^O{\'), ^~0(A2), !^^0(A4). (22) 

as ^ as as a^ 

We stress that this representation remains valid even in the case of large tan j3, up to tan (3 ~ 
0{mt/mb) - on the other hand, if the latter is not large, all the 1/1-4, wi and w^ can be 
dropped. Also, ai_s, xi, X2, yi and j/3 must be real since ra^guj-) are hermitian, while all the 
others can be complex. There are thus a priori nine new CP-violating phases entering in the 
squark soft-breaking terms, though only those of 04 and (less so) a^ and y^ are unsuppressed 
for moderate tan/?. 

Looking at the representation ( 12T1) . one may observe that, in contrast to the original 
formulation of Ref. [2], the symmetry principle which is behind MFV is not immediately 
visible anymore. Actually, the essential concept that the only source of flavour violation is 
the CKM matrix is still manifest in the Xj basis, and moreover: 

1. The symmetry principle which forms the basis of MFV is used once and for all, and 
leads to the basis shown, e.g. in Eq. (IT2l) . The relation between the standard MFV and 
the Xi basis can also be derived once and for all; it is given explicitly in Appendix [B] 
in terms of the running Yukawa couplings, to be defined in Sect. [31 

2. To make the most of the MFV principle, one still has to drop suppressed terms in 
the most general representation provided by MFV in Eq. ( IT5l) . As we have discussed 
above, it is easier and more transparent to do so in the Xi basis. 

3. The Xi coefficient^ have a direct, simple relation to the mass insertions [10], which are 
useful for phenomenological studies of flavour violations, as we will discuss in detail in 
the next section. 



4. The scalings of the Xi coefficients, as well as the decompositions ( l2T]) . are stable under 
the RGE's, as will be discussed in details in Sect. |3l Further, the formulation of the 
RGE's for the Xi coefficients is particularly convenient, especially because one can 
immediately see the order of the different terms appearing in their beta functions. 

2.4 MFV mass-insertions and their impact on phenomenology 

In many phenomenological applications, it is convenient not to assume anything about the 
structure of the soft SUSY breaking terms other than the established experimental fact that. 



"'in what follows we will often use the expression "the Xi coefRcients" or "the x^'s" to mean all Xi, yt and 
Wi coefficients. In order not to generate confusion, we will say explicitly when we mean specifically xi and 

X2. 



if present, they are almost diagonal. One then writes them as diagonal matrices plus small 
off-diagonal corrections and, when calculating observables, expands them in the off-diagonal 
terms - called mass insertions [TU]. Of course, in most practical applications, the exact 
diagonalization of the squark mass matrices is performed (see e.g. Refs. fTR IT2] for studies 
in the MSSM with MFV). But even if the use of mass insertions is thereby circumvented, they 
remain a very convenient tool to organize and identify possible sources of flavour violation. 
Indeed, the rich experimental information on flavour violations has been translated in bounds 
on the mass insertions [131 [E], and it is therefore useful to provide a relation between the 
MFV parameters and the latter. 

The LL and RR mass-insertions, defined as 

(X \IJ ^ ("^q) (xF ^/J ^ {mlY^ 

with F = U,D, are in terms of the MFV coefficients 





- V V* ""' + y*' 




(24) 
(25) 


'" *^ai|i/2ai + xi + yi + 2Rei/2^/' 
\ai 



up to completely negligible corrections of relative order A'^ or higher. In our approximation, 
all 6'^ji = 0, while only one 6^^ is non-zero, 

/rD \23 "^1 /rD \32* /0«^ 

i^RR) = 1^ 11/91^ , .. n/9 = i^RR) ' (26) 



^nf \ 10. 



but is nevertheless extremely suppressed since Wi ~ (9(A^°t|). The 5'f/(2)j^-breaking RL 



mass- insert ions are defined as 

^V 

(m2,)^^|i/2|(m|)'^-^|V2 ' ^"«^^ - I (mf,)^^ 1 1/2 1 (m^) ^-^11/2 



(^U N/J ^ ^"(^ ) (xD )IJ ^ '"dj-^ ) ^271 

\"RLJ I^™2 \/-7-|1/2I^™2 N.7jn/2 ' \"RLJ I /™2 N /-m /2 I ^™2 \ JJ 1 1 /2 ' ^ '^ 



and are 



s:(7 n32 _ T/ '^"'^4 _ Yj±( xU n31 



RL) 



[Orl) - VcbTZ- 



|a2|i/2|Si + xi + yi + 2 Reysl^/^ 



isi^r - v..^ ^""^ 



ai|i/2|a2|i/2 



■iD n32 _ T/ '^d'Ub _ Yj±(xD \Z1 



RL) 



^'^^^^''" |S3|V2|ai + xi + yi + 2Rey2P/2 ' ^^^^ 

while (5^^) ^2 _ ^^[7^)13 _ as well as (5^^)13 = (5^^) 21 = (5^^)12 = q under our approxima- 
tion of neglecting anything of order A^ or higher. The forms of the various mass-insertions 
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show that indeed, the CKM matrix elements still tune all the flavour transitions. In addition, 
the Xi basis permits to immediately judge of their respective strengths from Eq. fl22p . 

A prominent feature of Eqs. (^31 fIU\ [2H|) is the occurrence of several CP-violating phases, 
not related to the CKM one. Indeed, 1/2, uii, 0,4^, y^, W2, W4 can be complex, while all the other 
coefficients have to be real to satisfy the hermiticity of the squark mass terms. Though it is 
known that in the MSSM, MFV implies the presence of new CP-violating phases, up to now 
they have been considered only in the trilinear couplings. We see here that the parameter 1/2 
brings in an additional CP-violating phase in the LL sector alscO. Though it is competitive 
only for sufficiently large tan/3, since ?/2 ~ C'(A^t^), it could nevertheless play a role inb ^ s 
and b ^ d transitions, but not in s ^ d ones (which are in any case very constrained by Ek)- 
Since {SllY^ and {SllY^ are proportional to each other, it is not clear if, for example, this 
phase could explain the tension observed recently by the fit to the b ^ s transitions done in 
Ref . [T3] . Before drawing any conclusion about the validity of the MFV hypothesis at low 
energy, the role of these phases in the phenomenology should be investigated in detail. 

Concerning the CP-violating phases in the trilinear couplings, they can also have an 
impact on low energy observables, though mostly for b -^ s and b —>■ d transitions. Indeed, 
in the s ^ d sector, given the suppression of W2, even at large tan/5, the quadratic mass- 
insertions of the form (S'^l)^'^ (S'^l)^^* dominates |TH]. The CP-phase of 04 thus plays no 
role, and these transitions are entirely tuned by the CKM phase (this remains true when the 
contribution of the /x-term to the LR mass insertion is added, see Ref. [llj for more details). 
This is a priori not trivial looking back at the MFV parametrization ( IT5|) . because of the 
presence of the complex 6j and q terms in the expansion of A^. However, the contributions 
of these terms is always accompanied with light-quark masses, and is thus suppressed jTlj . 
This suppression is immediately visible in our parametrization, based on the counting rules 

m 



2.5 Lepton sector 

Since the RGB's for the (s)quark parameters depend also on the (s)lepton sector, we add 

m^ = ml [ae + bi^YlYe] , 
m| = ml [aj + 614 Ye Y^] , 
A^ = AoYe [as + feisYl Ye] . (29) 

For simplicity, we do not allow for flavour mixing in the lepton sector, and the counting 
rules describing the hierarchical lepton masses are also in powers of A as m,T-/mt ~ 0{X^), 
fTT'ii/'tTT't ~ ^{^^)y fTT'e/fnt ~ C^(^^)- Projecting the soft SUSY breaking matrices in the 
leptonic sector onto the Xi basis, we get 

m^ = Se + yeXi , m| = 07 + 2/7X1 , A^^ = asXi + W5X2 , (30) 



''The impact of having a complex 2/2 was partially analyzed in Ref. [17j . Indeed, in that work, though 
universality is imposed on the squark mass terms, so there is no new CP-phases, each trilinear term has 
a CP-phase at the GUT scale. As we will explore in some details later, the squark mass terms then can 
develop imaginary parts of precisely the MFV form by running down to the electroweak scale. 
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Relative to the leading terms, og (which incorporates Ur) as well as the j/j's can become of 
order one if tan/? ~ A~^: 

h^o{\%), y^^yi^oixHi), ^^o{x'). (31) 

04 Qq a-j ^ as 

Also in this sector, this representation remains valid in the case of large tan/5, up to tan/3 ~ 
0{mt/mb)) - on the other hand, if the latter is not large, j/ej can be dropped. 

3 Derivation of the renormalization group equations 

We have shown before that imposing the principle of MFV to the soft SUSY-breaking terms 
is, from a mathematical point of view, a mere reparametrization of these. One can therefore 
do an exact rewriting of the RGB's for the soft SUSY-breaking matrices into RGB's for 
the MFV coefficients, the Zj's appearing in Bq. flT^ . If the RG evolution does not make 
any of the Zi^s become huge, then we can safely drop also from the RGB's the irrelevant 
or redundant coefficients and get the RGB's for the reduced MFV parameters, the Xj's of 
Bq. (12T|) . In the following, we discuss this procedure in details, and give the RGB's for the 
reduced set of MFV parameters. The first question we have to address, however, is whether 
to allow our basis matrices (either the Y^^ and products thereof, or the Xi) to run or not. 
In order to do this, we first have to discuss how the Yukawa matrices (and correspondingly, 
the CKM matrix) run. 

3.1 RGE for the Yukawa matrices 

We first analyze the case of moderate tan/5, and discuss below the necessary modifications 
when tan/3 becomes large. Our starting point are the counting rules for the quark masses 
which we have defined in Bq. (1171) , the background values of the Yukawa matrices given in 
Bq. ([H]) and the Wolfenstein parametrization (I16p for the CKM matrix. We set the up-quark, 
down-quark and electron masses equal to zero and systematically neglect anything of order 
A^ or higher. The Yukawa couplings at the electroweak scale then have the following forms: 

Y„ (Mz) = VcXe + ytX, , Y^ (Mz) = ysX2 + y.X^ , Ye {Mz) = y^X^ + y.Xi . (32) 

For very large tan/3, these initial conditions may have to be amended, as discussed below. 

In the RGB for the Yukawa matrices themselves, but also of other matrices of the MSSM, 
products of several Yukawa matrices appear. But since the Xi basis is closed under matrix 
multiplication, all these RGB's corrections can be projected back on the X^ basis. Applied 
to the Yukawa matrices themselves, if we run according to the RGB's of the MSSM with the 
initial conditions in Bq. ( !32|) . we find out (as expected) that additional structures appear. 
But once our counting rules in A are enforced on the RGB alsqj, it turns out that it is 



^In view of the loop factor l/(167r^) in the definition of the beta functions, we keep leading terms only 
up to order A** in the beta functions, in contrast to A^ in the matrices. 
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sufficient to add only one term for Y^j and Y^^ to obtain an RGE invariant structure: 

Y„(/i) = ycil^)Xe + yt{iJ^)Xr^ + Ci(/i)Xi , 
Yd(/i) = ys{^)X2 + yb{^)Xi + Cb(/i)X5 , 

Ye(/i)=y^(/i)X2 + l/.(/i)Xi , (33) 

Tlie running of the three 3x3 Yukawa matrices collapses to that of only 8 independent 
parameters. 

We stress that the matrices Xj are held fixed - of course, the physical CKM matrix runs 
also, as one can easily realize by rediagonalizing Y^^^ (/^) in Eq. ( l33l) with q, cj, ^ 0. The 
CKM matrix at the electroweak scale is however given once and for all, and we use it to define 
a basis of numerical matrices, Eq. (1201) . in which to express the running Yukawa couplings. 
One can think of this basis as a fixed grid on which the RGE's for all fiavour-breaking 
parameters are projected. As will become clear in the following, this grid is particularly 
well-suited to enforce the MFV counting rules on the RGE, as it will permit to separate the 
rapid fiavour-blind evolutions from the much slower generation of fiavour-breaking effects 
through the running. 

The beta functions of the coefficients in Eq. flH^ . defined according to 



f4/3,. iV.16.^ (34) 

with t = \nQ/Qo, then read {yt = yt + Ct, yb = yb + Cb) 

Pyt = yt (Qut - ^«) + ChMt , Pet = Ct (Qyt - Ku) + ybVbyt , 

Pvi = yb {Qyl + yl- Kd) + ctym , Pc, = Cb {Qyl + yl- Kd) + Vbym , , . 

Py. = Vc (3y,2 - K.,) , /?,. = yr (%^ + ^yl -K,) , ^ ^ 

Pvs = Vs {Sg ^yl- Kd) , /?,, = y, {3y! + y'^ - K,) . 

Ku = ^gl + 3gl + ^^gl , Kd = K^- '^-gl , K, = Sgj + ^gl . (36) 

The coefficients Ct and Cb are zero by definition at /x = Mz, but are generated by the running 
at any other scale (see Fig. [5] in Appendix El)- The leading terms in their beta- functions 
are yfyt and y^yb, respectively, and from these we can infer (albeit with some degree of 
arbitrariness) how to count them. In the following we adopt as counting rule: 

Ctr^OiylX^) and c,r^O{y,\') . (37) 



where 



According to this, the particularly simple beta-functions in Eq. flHoj) are accurate up to 
corrections of order A^ (and in many cases to better than this). We checked numerically 
that indeed, Eq. ( !35|) agrees with the full one-loop running of the Yukawa couplings, after 
projecting them back on the basis of Eq. ( I33l) . to better than 0.1% for moderate tan/3. 

At large tan/3, the Yukawa couplings at the electroweak scale are not always simply 
related to the quark masses, and the initial conditions at /i = Mz in Eq. (!32l) have to be 

12 



corrected. Indeed, non-holomorphic Higgs couplings arise beyond leading order from the 
combined breakings of the U{1)pq symmetry by the // and b terms, and of supersymmetry 
itself by the soft-breaking terms [20]. After electroweak symmetry breaking, non- flavour 
diagonal, tan /5-enhanced contributions to the down-quark mass matrix emerge. The net 
effect is that in the basis in which Y„ is given in terms of the physical up-quark masses and 
CKM matrix, Y^ is not diagonal: 

Yu = X^V , Yd = Ad(l-A)-' , (38) 

where VuXu, VdXd are the (diagonal) physical quark mass matrices, and V the physical CKM 
matrix^. All the loop-induced, tan /9-enhanced corrections are in A, which has the MFV 
expansion (under the assumption that Y^ is still sufficiently hierarchical) 

A = tan/3 (sqI + ^lY^Y^ + saY^Y, + ^sY^YrfYtY^ + e^YlY^YlY^) . (39) 

The Ei parameters are loop suppressed and depend on the mass-spectrum of the specific 
MSSM model under consideration. However, for large tan/?, this suppression can be com- 
pensated, at least partially, leading to a relation between Y^ and A^ with corrections which 
may become of 0{1) if ejtan/? ~ 1. 

To see what could happen in that case, we remark first that the Cayley-Hamilton identity 
([9]) implies 

6X2-6(X)X + 3((X)'-(X2)) 

2 (X3) - 3 (X) (X2) + (X)= 

For X = 1 — A, this becomes 



x-^ = ...J .,...:!/. 's^ ■ (40) 



A)-' = 7],l + 7]2A + r]3A\ (41) 



with the coefficients 



V^/V3 = 1 - (A) + i ((A)^ - (A^)) , r^^M = 1 - (A) , (42) 

%' = 1 - (A) + ^ ((A)^ - (A^» + I (A^> (A) - I (A^> - I {Af . 

There is no approximation in this formula: Cayley-Hamilton allows to completely resum the 
geometric series expansion of (1 — A)~ . When the mass spectrum is such that ejtan/3 is 
small, these formula can be expanded to first order, leading to Y^ ?« A^ (1 + A). 

Provided the £« are not too large, the coefficients rji are 0{1) numbers. In that case, the 
expansion (1391) is certainly valid, and all the tan /3-enhanced corrections can be absorbed as 
shifts in the values of i/b (Mz), Vs (Mz) and c^ (M^), since (1 — A)~ has the same form as 
the RGE effects. We stress that, except for very large St, the counting rules are not upset 
by these shifts. Only Cb tends to become somewhat larger than before (see Appendix |A] for 



^Tliis quark basis is different from the one of e.g. Ref 2^, where the background value for Y^ is kept 
diagonal, but at the cost of having Yu = XuV with V' different from the CKM matrix. Here, the down 
quark fields are already mass-eigenstates, since once loop corrections are added, M^ — VdY a (1 — ^) = v^Xd- 
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a numerical application), but this increase is quite mild, and the RGE's, Eq. (1351) . remain 
accurate to better than 1%. 

In the present paper, our aim is to probe the behavior of the RGE's also in the large tan/? 
scenario. Since the Si depend on the mass-spectrum, which itself depends on the matching 
conditions at the electroweak scale, and hence on the £j, this problem can be properly solved 
only by setting up an iterative procedure, and this is beyond our scope. Our analysis could 
in fact be viewed as the first of these iterations, allowing to derive an approximate mass 
spectrum, and to compute the £j parameters. Even if in the end the £« will turn out to be 
very large, such that r^j S> 1, we stress that the counting rule method developed here and in 
the previous section could be generalized to more complicated initial conditions, simply by 
allowing for additional structures (and parameters) in the Yukawa couplings in Eq. (l32l) . 

3.2 RGE for the MFV parameters 

We now consider the soft SUSY breaking terms. The beta functions for them are known 
(even up to two loops, [18j) and can also be projected onto the Xi basis. In principle, one 
could generate new structures beyond those given in Eq. fl2T|) if the beta functions of some 
of them were of a different order than the coefficients themselves. We have verified that this 
does not happen, and that the structure in Eq. fl^ is indeed RGE invariant if we stick to 
our counting rules {i.e. deviations from this structure arise in the running from contributions 
to the beta functions which are suppressed by at least A'^). We stress that this is only a 
necessary (but not sufficient) condition to ensure that the MFV hypothesis is RGE invariant 
- an exponential growth can of course be generated by a beta function which is of the same 
order as the parameter itself. We will then have to study the behaviour of the solutions of 
the RGE in order to reach the desired conclusion. 

The beta functions of the new coefficients turn out to be remarkably simple. They read 
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(43) 



for the squark soft-breaking terms: 

/3^i = 2|/2 (m%^ + ai + a2 + xi + X2 + Ret/s) + 2 (|a4|2 + l^gp) , 
Pyi = 2|/fe ("^L + 5i + as + yi + Re?/2 + Vs) + 2 (iSsp + \y^f) , 
/^y2 = ^6(3^1 + ^2) + y! (l/i + 1/2) + 2 (a^i/5 + a^yl) , 

vr.2 . /3 _ 32^2|/\/f 12 32^2|/\/f 12 4^2 

P^2 = 4|/2 (77^2^^ + a^ + ^2 + a;i + a;2 + t/i + 2 Ret/2) + 4|a4 + |/4p , 

m2 ■ /^ — 32 2|/\/f 12 8 ^21 /\^ 12 , 2 2 

/5y3 = 4|/^ (m|^^ + ai + d3 + xi + yi + 2Rey2 + ys) + 4|a5 + 1/5P , 
(3^, = 2wiyl + 4 (55 + 1/5) wl - AA\^ ((xi + 1/2) y^l/s + W3I/5) , 

A^ : 13a, = ^4 (181/2 - /CJ + y, {lly^y^ + 2|/5|/fe + iT^ , 

/5y4 = 1/4 (l/fe + ^y'i - Ku) + a^yl + 2a5|/b|/t , 

/5«,2 = ^2 i^yt - Ku) + yc (Qyt (^4 + y^) + K'J , 

A^ : f3a, = S5 (18t/2 + t/2 _ J^^) + ^^ (11^^^^ + 2y4t/t + 2fi8yr + K) , 

/3y5 = 1/5 (7l/b + 1/i^ + y? - ^d) + a^y^ + ^a^ym + c^K'^ , 
Pw3 = W3 {3yl + yl- Kd) + ys (Qyt (05 + I/5) + '^agyr + K'^) , 
/3^4 = W4 {8yl + y'^ + yl- Kd) - AX'^ {2y,yt (04 + t/4) + wsl/?) , 

and for the slepton soft-breaking terms: 

mi: f3d, = -Qgimi' - Igfmi' - IgfS , 

j3y^ = 2yl {m]j^ + ae + 07 + ye + yi) + 2|a8p , 
m|: /5,7 = -f^?|Mip + |^?5, 

I3y, = Ayl {m]j^ + ae + 07 + ye + yi) + 4|a8p , 
A^ : 13d, = ds {y2yl + 'iyl - K,) + y, (6 (fig + a:io) yh + i^e) , 

/5«;5 = ^5 (y^ + 3^^^ - i^e) + 1//. (6 (as + 1/5) yb + 2~a^y^ + iT^) , 

where Vch — AX^ ~ 4-10^2 is the element of the CKM matrix, and where we have introduced 
the abbreviations 

K = f ^1^3 + 6^2'M2 + ^alM^, K'd = K + ^gfM^, K = Q9IM2 + y^?Mi , (45) 
S = mjj^ - mjj^ + 3 (ai - 2a2 + 03 - Oe + 07) + a;i - 2^2 + yi + 2 Re 1/2 + I/3 - ye + I/7 • 
The only other RGE's depending on the sfermion soft-breaking terms are the Higgs param- 



(44) 
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eters, which now take the form 

/3^ = /i {3y^ + 3yl + y^ + 3yl + yl- 3gl - ^l) , 
A = b {3y^ + 3yl + yl + 3yl + yj - 3gl - Igl) 

+ /i (6 (04 + 1/4) yt + Q (as + 2/5) yi, + 2asyr + Qw^ys + Su^sy^ + 65-2^2 + IglMi) 
j3^2 = Qyl (mjj + fii + 02 + a;i + X2 + z/i + 2 Re ^2) + 6|a4 + ^4^ 

-Hit \ ti / 

- Qgl\M2\^ - -gl\Mi\^ + -g^S + Qy^, {m\^ + Si + ^2) + 6|i/;2|' , 
/?m|, = 62/6 ("^Hd + fli + 03 + xi + yi + 2 Re?/2 + Z/3) + 6|a5 + ^5^ 



'd 



- 6(?2'|M2|' - -gi\M^\^ - -g\S + 2y2 {rn\^ + 5^ + 07 + 1/6 + 7/7) + 2|58r 

+ 6?/2 (ml,^ + ai + 03) + 6|w;3p + 2^^ (m^,^ + fig + ^7) + 2|w;5|' • (46) 

In these last equations, corrections of relative 0[\^) were kept. They can be relevant at the 
percent level in some corner of parameter space because of their rather fast running. 

We stress that the RGE's (H3l HH HE]) are to be taken as they are in the case of large 
tan/5, i.e. when t^ ~ (9(A~^), as long as the Yukawa couplings have the structure shown in 
Eq. (l33ll . On the other hand, for tan/3 of order A~^ or even A~^, some couplings and the 
corresponding beta functions change their order and may become negligibly small. In that 
case, the flavour structure of the theory as well as the RGE's become significantly simpler. 

4 Running MFV in the moderate tan/3 case 

In this and the following section, we will illustrate the behaviour of the various parameters as 
functions of the scale with the help of a few examples. We first discuss the case of moderate 
tan /3, and then the case of a large one. In order to do this, we take as reference two of the 
Snowmass benchmark points ^19], and use the MFV parameters to explore flavour violations 
around these points. These numerical examples are essential in order to understand the 
solutions of the RGE's. In case of moderate tan/3, however, the structure of the MFV 
MSSM, and of the corresponding RGE's, simplifies so much that one can provide a semi- 
analytical solution of the RGE's. We will now first derive these, then present the numerics, 
and show how one can understand the observed behaviour with the help of the analytical 
formulae. 

4.1 Analytical solutions 

4.1.1 Solutions for the a^'s 

We first observe that the term proportional to S in the RGE's for 01,2,3 is typically very small: 
it is multiplied by the small gauge coupling g\ and a small coefficient, and the combination 
of massive coupling constants which appears in there is zero at the GUT scale for initial 
conditions of the mSUGRA type. The numerical examples we will discuss below will show 
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that this term can indeed be safely neglected, unless one takes very special initial conditions. 
Dropping S, the RGB's for 01,2,3 admit a simple explicit analytical solution (to = ^(Mgut)): 



5i(t) = ai(to) + - (|M3(t)r - |M3(to)r) - - {mit)\' - \M,{to)n 
(|Mi(t)p-|Mi(to)P) , 



198 



a^it) = a,{to) + - (|M3(t)P - |M3(to)r) - ^ (|Mi(t)p - |Mi(to)r) 
asW = dsito) + I (|M3(t)r - |M3(to)r) - ^ (|Mi(t)p - m{to)\') 



(47) 



Typically, the term |M3(t)p, which grows fast when evolving down to the electroweak scale, 
will dominate over the rest, so that all three Sj's turn out to be of the order of the gluino 
mass squared at the electroweak scale. 

The RGB's for 04,5 simplify as follows in the moderate tan/3 caso 



Pa, = S5 {Wb + yl- Kd) + VhK'^ . 



Both RGB's are of the form 



dA 



fA{t) + A{t)gA{t) 



dt N 
with /^(t) and gA{t) known functions. The solution of this equation reads 

"to 



A{t) 



Ait.)-^l 



t GA{t') 



GA{t) , GA{t) = exp 



1 r° 
— / dt'gAit') 



(48) 



(49) 



(50) 



4.1.2 Solutions for m?/ and xi 2 

The RGB's for m'jj^, xi and X2 are coupled, even in the moderate tan/3 case. In that case, 
however, one can diagonalize the beta functions and bring them to the form ( H9l) . by taking 
the following combinations: 



rii 



[mjj^ + xi+ X2] 



^2 



2 L "" ^ 2 

The / and g functions for these combinations read 



[m^^ -xi- X2] , ns = 2xi - X2 ■ 



(51) 



^ni = 12l/t 

9n2 = 
9n, = 



fn, = -Sgl\M2\^ - -gl\Mi\^ + 6 [y^{di + fia) + 1541'] 



Jn2 



-H\M2?-^-gl\M,\' 







*Note that the imaginary parts of 04,5 miss the term proportional to K'^ ^ in the beta function. 



(52) 
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and are all known functions. One can plug them into the general solution (!50l) and gets: 



ni{t) 



ni{to) - — at 



NJt G^,(t') 



GnM 



n2{t) = n,{to) - I {m{t)\' - |M2(to)r) - ^ (|Mi(t)p - |Mi(to)r) , 

nsit) = ns{to) , (53) 

and from these solutions, reconstructs the physical parameters: 

^Hu =ni+n2 ^1=0 [^1 -^2 + ^3] X2 = - [2(ni - ^2) - na] . (54) 

The most important feature of these analytical solutions is that the integral containing the 
function /„^ (which gives a negative contribution) dominates, because it contains the fast 
growing functions Oj, which behave like the gluino masses. This growth is only partially 
compensated by the decreasing function 0^ (as t decreases from to down to the electroweak 
scale). The function n2, which is positive, grows much more slowly, approximately linearly 
in to ~ t. From these analytical solutions, one can also clearly see how much the values of 
these parameters at any scale depend on their initial conditions at the GUT scale: 

ml^t) = ^^i^ [G„,(t) + 1] + 1 (xi(to) + X2{to)) [G^At) - 1] + . . . 

XI (t) = ^^%^ [G„,(t) - 1] + ^ [G„,(t) + 5] + :^ [G„,(t) - 1] + . . . 
bob 

X2{t) = ^^%^ [G„,(t) -l] + ^-lM [G„^(t) _ 1] + :^ [G„^(t) + 2] + . . . (55) 

where terms which are independent from the initial conditions have been omitted. In order 
to understand this dependence, one only needs to know the behaviour of the function Gm ■ 
As just mentioned, this decreases from the GUT to the electroweak scale, and in fact does 
so almost linearly from 1 to about 0.23. 

For moderate tan/3, the RGE for mfj can be solved in terms of known functions, whose 
evolutions are not influenced by mj^ itself. The / and g functions read: 



9ml = ^yb + '^yl ^ (56) 



4.1.3 Solutions for the yj's and the Wj's 

The same applies to all the remaining parameters. One can give their solution in the form 
of Eq. (l50l) simply by identifying the / and g functions for each of them. We give here a 
couple of examples: 

9yi = 2l/b , fy, = 2y^{mjj^ + Si + as) + 2|S5|2 , 
9y2 = Vt > fy2 = Vb^i + VtVi + 2(042/4 + «52/5) • 



In all these cases, the analytical solutions allow one to immediately understand what depends 
on what, and in particular, which initial conditions influence the behaviour of the various 
parameters at the electroweak scale, as we are now going to show in the following subsection. 

4.2 Numerical example: the SPS-la benchmark point 

The SPS-la point is specified by tjiq = —Aq = 100 GeV and mi/2 = 250 GeV at the GUT 
scale, and tan f3 = 10. The running of the MFV parameters in the neighbourhood of the SPS- 
la benchmark point is illustrated in the Figs. [H We have checked that the running evaluated 
by solving our simplified RGB's gives almost identical results to the running evaluated with 
the full one-loop RGB's of Ref. [18] (for the numerical analysis and the figures, we have 
included also the A^ corrections to the beta functions given in Appendix [C]) . We have solved 
both sets of differential equations with the same input at the GUT scale and then compared 
the full mass matrices at the electroweak scale - the differences in all the entries are always 
consistent with our counting rules, i.e. of the order of the terms neglected in our simplified 
RGB's. 

The most prominent feature emerging from the figures is that all the parameters ratios 
tend to converge to one point at the low scale, independently of the starting point at the high 
scale - i.e. they show a fixed-point kind of behaviour, as recently observed in Ref. [H]. We 
stress, however, that by separating the leading parameters (the Oj's and the Xj's) from the 
suppressed ones, we can better see the behaviour of the small flavour violations, governed 
by the ?/j's and the Wj's. Our analysis shows that these also converge to fixed points. 

The analytical solutions discussed above allow us to clearly understand the origin of these 
fixed points. First of all, we stress that indeed the running of the Oj's is almost independent 
from the Xj's, (which come in only through the combination S), and is strongly dominated 
by the term proportional to the gluino mass and strong gauge coupling, cf. Bq. (H71) . Table [1] 
shows that the Sj's tend to the value 8|M3p/9 (for = 1,2,3) at the low scale, up to small 
corrections, which of course can be calculated exactly with Bq. fH7|) . For 04 5, we do not have 
an explicit analytical solution, but from Bq. fl^Ujl we see that the integral over the function 
/ tends to a value of the order yt^Msit) at the low scale. For 04 (05), the integral Ga^^^^ 
decreases (increases) from 1 to about 0.4 (3.3) in going from the GUT to the electroweak 
scale, and so we understand why the proportionality factor is smaller (larger) than one at 
the electroweak scale. 

The analytical solutions for the Xj's, Bqs. (153113^ . also explain the corresponding numbers 
in Table [U and Fig. [H First of all, the values of xi and X2 at the electroweak scale, with 
initial conditions a;i(to) = 2^2(^0) = 0, satisfy the relation 2x1 (t) = 2:2 (^) according to the 
solutions ( l53l [5^ - a comparison to the numerical solution of the exact RGB's shows that 
this is indeed verified to better than one per mijj. We can therefore consider the solution of 
only one of the two, which reads: 



m, 



2 



GnM f\jnAt') - fn,{t') , 1 Z"*" 



Xl(t) = ^ [G„,(t) - 1]-^^^/ dt' ^^^^'-J, /T^ +irr dt' fn.it') 



6 ^ "^^^ ' 3N J, G^,{t') 3NJ, 



Gni{t) 
Gni{t')_ 



■ (58) 



'Notice that a;2/a2 = Ixxja^ = 2x\la\ ■ ai/a2 = —0.36 ■ 1.06 = —0.38. 
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Figure 1: Upper four plots: RGE evolution of the MFV parameters for the SPS-la benchmark 
point. The solid curves always show the evolution with mSUGRA type of initial conditions. 
In the upper-left panel, the dashed (dotted) lines show the evolution when a;i(2)(MGUT) = 
±Mq. In all other cases, for each parameter, only three curves are shown - the upper 
and lower ones (always shown as dashed) correspond to different initial conditions for that 
single parameter. Lower two plots: the real and imaginary parts of the mass insertions b\(2) = 

i^RL )^^/^s = i^RL j^^/^id- For these plots, the initial conditions are varied independently 
for A*^ and A-^, allowing for a CP-phase as Aq = re**^ with r = —>■ 200 GeV and between 
±180°. Dashed lines in the 62 plot show the impact of the initial conditions for A*^, while the 
sensitivity of 61 to those for A^ is negligible. The behaviours of the other mass- insert ions, 
Eqs. (I2S1 [2ni [2H]) , can easily be obtained from those of the parameters shown in the upper four 
plots. In particular, note that j/2 (Mz) <^ Xi (Mz), both because of the tan/5 suppression, 
and because of RGE effects. 
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Mgut 


Mz 


ai/M| 

52/M| 

as/Ml 

aJ{-ytM^) 

a^/i^-ybM'i) 


0.16(1 + Ai) 
0.16(1 + A2) 
0.16(1 + A3) 
0.40(1 + A4) 
0.40(1 + A5) 


0.86 + 0.019A1 
0.8I + 0.015A2 
0.8O + 0.019A3 
0.8I + 0.03A4 
1.48 + 0.14A5 


xi/ai 
yi/ai 
y2/ai 


'^1 
e2t}X' 


-0.18 - O.OO3A2 - O.OIA4 + 0.025i - 0.003(52 
(-7.6 + O.lAi - O.IA3 - I.IA5 + 0.3ei) ■ lO'^ 
(4.4 - O.lAi + I.IA4 + I.4A5 - 0.45i - 0.4ei 
+3.0e2 - 2.5e4 - 4.8^5) ■ 10"^ 


X2/a2 


S2 


-0.38 - O.OlAi - O.O2A4 - 0.01^1 + 0.02(^2 


y^/as 





(-1.53 - O.O2A1 + O.OlAi + O.O2A3 + O.O3A4 - O.2A5 

-0.025i + 0.04e3 - 0.2^5) ■ 10-2 
~ -10-6 


2/4/04 
102/04 


ctlyt + ^At^X^ 
yjyt + ^2 A'' 


(-4.1 - O.IA4 - O.6A5 + 2.3e4) ■ 10-3 
(5.1 -O.4A4 + 0.7772) -10-3 


1/5/05 
ws/as 
104/0,5 


Cfe/z/b + S5 

ys/yb + ?73A^ 


-0.079 - O.OI2A4 - O.OO3A5 + 0.084^5 

0.019 - O.OO2A5 + 0.005r/3 
(1.0 + O.IA4 - O.IA5 + 0.2r/3 + 2.3774) ■ 10-^ 



Table 1: Dependence of the MFV parameters on the initial conditions when these are taken 
in the neighbourhood of the SPS-la benchmark point. Any correction below one percent 
has been omitted. 



At the electroweak scale tew the three terms contribute, respectively 

a;i(tcw) = -(1.3 + 73.1 + 2.9) ■ 10^ = -77.2 ■ 10^ GeV^ 



(59) 



which shows that the dominating contribution comes from the integral over /„j — /„2 which 
contains only the Oj's. In the ratio x\/a\^ this is the part which tends to a small constant 
value at the electroweak scale, which almost looks like a fixed point. To understand why 
this almost looks like a fixed point, we have to analyze the dependence of a;i(t) on the 
initial conditions. The relevant formulae have been given above, cf. Eq. (13^ . and show 
a linear dependence, with coefficients of order one. If we set, e.g. X\{t^ = di{to)6i and 
3^2(^0) = d2(to)62, as indicated in Table[Il then at the electroweak scale the effect is equal to: 



Axi 

AX2 
0.2 



ai(to) G„^(tcw)+5 02(^0) Gni{tcw) ~ ^ 



ai(te 



6 



Si^ 



ai(to) Gni{tc 



1 



Ol{tcw) 
0-2 (to) Gniitc 



6 



(60) 



02 (te 



'2- 



0.2(4 



where Axi = Xi{te-„)\s ^p — Xi{te„)\g ^^, and inserting the numbers Gni(4w) = 0.228, and 
oi(to)/ai(tew) = 0.023, 02 (to) /a2 (tew) = 0.024, vahd for the initial conditions specified in 
Table [H one perfectly reproduces the numbers there. Similarly one can explain why the 
sensitivity to the initial conditions of the a^'s is small: these enter through the second 
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integral in Eq. (l58l) and contribute as follows. If we take aj(to) = ^^oi^ + \) for ^ = 1, 2, 3 

and a4(5)(to) = yt(b){to)A{^ + ^4(5)) then we have (for Axi = a;i(tew)|Afc^o - a;j(tew)|A,,=o) 

=-m^[0.13(Ai + A2) + 0.38A4] . (61) 

Dividing these numbers by ai(tew) one reproduces the numbers in Table [H Note that the 
sensitivity of xi/fii to Ai receives another contribution due to the denominator dependence 
on Ai. This is easy to calculate and brings the total coefficient below 1% - for this reason a 
term with Ai does not appear in the line for Xi/ai. We stress again that while the sensitivity 
of Xi itself to the initial conditions of the parameters which appear in its beta function is 
described by a coefficient smaller than one, but not tiny, the ratio xi/ai has tiny coefficients 
because ai(to)/Si(tew) is tiny for the initial conditions given by the SPS-la benchmark point. 

All other coefficients appearing in Table [1], which specify the sensitivity of the various 
parameters to the initial conditions of all the others, can be understood analogously on the 
basis of our approximate analytical formulae. 

Finally, concerning the CP- violating phases, we see from Fig. [1] that if 5^^ or S^j^ have 
an imaginary part at the GUT scale, even if large, it is much suppressed by the RGE effects. 
This is easily understood on the basis of the simplified analytical solutions: the imaginary 
part of 04 misses the piece of the beta function which is responsible for the fast running of 
the real part, cf. 0481) . At the electroweak scale, only a residual phase remains, which acts as 
a small correction to the CKM one present in Vts or Vtd- Note that at moderate tan/3, (5^^ 
is significantly affected by the initial conditions for 5^^: a phase for the latter can generate 
one for the former. On the other hand, any phase present in the trilinear terms has only a 
very limited effect on 1/2, relevant for {6ll)32 and {Sll)3i mass- insert ions. This follows also 
from Table [T], where it is apparent that the initial conditions for the trilinear terms affect 
only mildly the final value of 1/2 at the electroweak scale. Conversely, if 7/2 (to) has a large 
imaginary part, its running is very similar to the real part, i.e., it also decreases dramatically 
(see upper right plot in Fig. [1]). Further, looking at the RGB's, ?/2 cannot communicate its 
phase to the trilinear terms. Therefore, all in all, if the SPS-la is extended to include as 
many CP-violating phases as allowed by MFV at the GUT scale, all these phases are strongly 
suppressed at the electroweak scale. In that case, MFV essentially collapses onto the real 
parameter case, usually considered in the literature (see e.g. [I2]), where all CP-violating 
observables are tuned entirely by the CKM phase. 

4.3 Running from the electroweak up to the GUT scale 

The numerical analysis of the SPS-la benchmark point has shown that if one assumes that 
the MFV hypothesis is valid at the GUT scale, one has rather definite predictions at the low 
scale. We can now turn the question around and ask how MFV evolves towards the high 
scale if one assumes to know the parameters at the low scale. In particular, it is interesting to 
analyze what happens if the boundary conditions at the low scale are chosen such that they 
are rather far from the "fixed point" discussed above, but still well within the naturalness 
assumption of MFV. 
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Figure 2: RGE evolution of the MFV parameters for the SPS-la benchmark point, if the 
boundary conditions are fixed at the low scale. The dashed (dot-dashed) lines correspond 



to the change of boundary conditions Xi(^2)iMz) = —x' 
plot reach Xi/ai ~ 18 and 0:2/02 ~ 40 at the GUT scale. 



SPS-la 
1(2) 



{Mz). The curves quitting the 



We illustrate this in Fig. [21 in which we use as boundary condition at /i = Mz all the 
values obtained from the evolution of the SPS-la point, with the exception of Xi 2 for which 



we use xi^2{Mz) 



-X 



SPS-la 

1.2 



{Mz). The figure shows that a point far from the "fixed point" 



at the low scale evolves to very high values of the MFV parameters at the high scale. This 
behaviour is easily understood on the basis of our analytical solution: from Table [1] we read 
off that a change of natural size of the initial conditions for xi produces a change of about 
15% at the low scale. So, if we want to induce a change of 200% at the low scale, we need to 
make a change of more than one order of magnitude at the high scale. One could imagine 
changing more parameters at the same time in such a way that the evolution would be of 
MFV type all the way to the GUT scale. Table [1] shows, however, that this is not possible. 
Changing any of the other parameter around the SPS-la benchmark points, while staying 
in a range compatible with MFV, does never generate a large change for Xi/ai at the low 
scale (similar conclusions can be reached for most other parameters). 

The analysis of the previous section points to a possible way out: if the ratio ai(to)/Si(tew) 
were of order one instead of 10~^, much larger changes in the low energy MFV parameters 
would become possible. This ratio, however, is completely fixed by the RGB's and initial 
conditions, cf. Eq. (H71) and can be changed significantly only by changing ai(to) with respect 
to M3(to) such that the former becomes at least as large as the term proportional to the 
gluino mass squared in Eq. (H7I) . This, however, brings us far away from the SPS-la point, 
and moreover, will tend to make the squarks heavier, which makes their contribution to 
low-energy flavour violations, and correspondingly the phenomeno logical interest, smaller. 
Such a situation is similar to the one realized in the case of the SPS-4 points (independently 
of the size of tan/5), as we will discuss in the next section. 

We conclude that, if the spectrum at the low scale corresponds to the SPS-la input at 
the high scale (note that the RGE's discussed here show that an MFV-compatible change of 
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the boundary conditions at the high scale has barely any influence on the spectrum), and if 
the values of the MFV parameters measured at the low scale are far away from those given in 
Table [H then the MFV hypothesis has to break down somewhere before reaching the GUT 
scale. 

This may indicate that either MFV has emerged from the RGE evolution of the pa- 
rameters, starting from a non MFV kind of MSSM, or that the MFV parameters have a 
dynamical origin at a scale much lower than the GUT scale. The former solution cannot 
be really investigated with our simplified RGB's, because if any of the MFV parameters is 
far from its assumed order, our RGB's lose their accuracy. One would have to study this 
solution with the full RGB's, and see how much fine-tuned it is, or whether it can viewed as 
natural in any sense. We leave this question open. 

5 Running MFV in the large tan/3 case 

5.1 Analytical solutions 

When tan/? ~ 1/A^, the parameters Hi become in principle of order one and cannot be ne- 
glected. The system of the RGB's is then not amenable anymore to an analytical solution. 
The differential equations are all coupled and a diagonalization, although possible in prin- 
ciple, is too complicated to be of any use. Therefore, we start directly with the numerical 
analysis and make some remarks about how one can understand Table [2] on the basis of the 
analytical solutions provided in the previous section, plus some corrections specific to the 
large tan/3 case. 

5.2 Numerical example: the SPS-4 benchmark point 

The SPS-4 benchmark point has ttlq = 400 GeV, mi/2 = 300 GeV, Aq = aX the GUT 
scale, and tan/? = 50. With these parameters, too large supersymmetric contributions to 
b —>■ s'j are generated, such that the total rate comes out much lower than the measured 
one - this point is therefore not compatible with the phenomenology at the low scale. We 
use it nonetheless for illustration purposes. Another important aspect is that Aq is set to 
zero at the GUT scale, and so a^ = d^ = there. If we want to vary the flavour violating 
parameters in the trilinear couplings, we flrst need to deflne a scale of naturalness for them. 
For this we use the value of d^ and as at the low scale and set at /i = Mqut 

yi^-ytM^tl\\ W2^-ytM^X'' , (62) 

and 

1/5 ~ -i/feMg , W3 VbM^X^ , W4 ~ -ybMsX"^ . (63) 

As soon as we make any of these parameters different from zero, their ratio to 0,4 or 0,5 at the 
high scale therefore diverges, as it is seen in Fig. El Running down to the low scale, however, 
the ratio quickly converges to its natural range, and becomes one or smaller. 

The fact that tan/5 is so large means that most of the yi and Wi parameters change 
their order, since A^tan/? ~ 0.6. The curves in Figure [3] show that the behaviour is indeed 
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Figure 3: Upper four plots: RGE evolution of the MFV parameters for the SPS-4 benchmark 
point. The solid curves always show the evolution of the mSUGRA type of initial conditions. 
In the upper-left panel the dashed (dotted) lines show the evolution when a;i(2)(MGUT) = 
±Mq. In all other cases, for each parameter only three curves are shown - the upper and 
lower ones (always shown as dashed) correspond to different initial conditions for that single 
parameter. Lower two plots: the mass- insert ions 5i(2) = {^^j^ )'^^/Vis = {^rl )'^^/^rf) with 
the initial conditions at the GUT scales varied as explained in the text, but allowing in 
addition for a large CP-phase (between ±180°). In this case, b\ ~ 04 is entirely radiatively 
generated, since 04 = is set to zero at the GUT scale. The behaviours of the other mass- 
insertions, Eqs. fl2^ [2Sl [2S|) . can easily be obtained from those of the parameters shown in 
the upper four plots. In particular, note that though ?/2 iMz) is still smaller than x\ (Mz), 
it is much less suppressed by RGE effects than in the SPS-la case, and so is its CP-violating 
phase. 
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Mgut 


Mz 


ai/M| 

52/M| 

as/Ml 

aJ{-ytM^) 

d5/{-ybM3) 


1.8(1 + Ai) 
1.8(1 + A2) 
1.8(1 + A3) 

A4 
As 


1.06 + O.2IA1 + O.OIA2 + O.OO6A3 

1.01 + O.O2A1 + O.I7A2 + O.O2A3 + 0.008((5i - 2S2) 

1.00 - O.OlAi + O.O2A2 + O.2IA3 - 0.004((5i - 2S2) 

0.78 + O.O9A4 + O.O2A5 - 0.09e4 - O.OS^g 

0.93 + O.O2A4 + O.I6A5 - 0.03e4 - 0.17^5 


xi/ai 
yi/a-i 
y2/di 


<5i 
e2t^A6 


-0.21 + O.OlAi - O.O2A2 - O.O2A4 + 0.185i 
-0.02^2 + O.Olei + 0.02e4 - O.Ol^s 
-0.16 + O.OlAi - O.O2A3 - O.O3A5 + O.OWi 
+0.07ei - 0.01e3 + 0.025^ 
(0.6 - O.lAi + O.6A4 + O.7A5 - lA5i - 0.9ei 
+6.3e2 - 0.9e4 - 1.6^5) ■ lO'^ 


X2/a2 


S2 


-0.43 - 0.04Ai + O.OIA2 + O.OIA3 - O.O3A4 + O.OIA5 
-0.05^1 + 0.15(52 - 0.02ei - 0.04e2 - 0.01e4 + O.Ol^s 


1/3/03 
wi/a3 


estlX' 



-0.32 - 0.046Ai + O.OI2A2 + O.O2A3 + O.OIA4 - O.O4A5 
-OM61 + 0.01(52 - 0.02ei - 0.03e2 + 0.06e3 - OM65 
-10-5 


1/4/04 
^2/54 


e4t2A6 (*) 
V2X' 


-0.09 - O.O2A4 - O.O4A5 + 0.14e4 + 0.02(55 

(5.0 - I.IA4 - O.IA5 + 0.364 + 1.9r]2 + 0.3^5) ■ 10-3 


1/5/05 

Wg/ag 

^4/05 


VsX' 

V4X' 


-0.10 - O.O46A4 - O.O23A5 + 0.015e4 + 0.36(55 
0.023 - O.OO5A5 + 0.001e4 + 0.003^5 + 0.026r/3 
(1.2 + O.3A4 - O.2A5 + 0.1e4 + 0.2^5 + 1.3r73 + 7.8774) ■ 10"^ 



Table 2: Dependence of the MFV parameters on the initial conditions for the SPS-4 point. 

Any correction below one percent has been omitted. 

(*) Since 04(5) (Mgut) = 0, at this scale we normalize 7/4 and W2 (1/5, W3 and W4) by —ytM^ 

(-Z/6M3). 

different than for the SPS-la benchmark point. Since now more parameters are of order one 
(or almost), the mutual influences of the various parameters are more important, though 
still quite small, as can be seen in Table O The most prominent difference, however, is that 
the "fixed points" are now much weaker - the value of the parameters at the low scale is 
influenced more sensitively by the starting values at the high scale. The same two statements 
can be made for the imaginary parts of the complex coefficients: their suppression through 
RGE effects is now much weaker, and the presence of a CP-phase for one parameter is felt 
more strongly by the others. For example, if A*^ and/or A^ involve complex phases at the 
GUT scale, 2/2 does develop a small but significant phase at the low scale. 

Can we understand why this happens? Also, and more importantly, is it the large tan/? 
or the large mo which is responsible for the different behaviour of the SPS-4 benchmark 
point? A first clue is given in Fig. |H where the boundary conditions are set according to 
the SPS-la, but for tan/3 = 50. Obviously, the fixed point behaviour is to a large extent 
unaffected by tan/?, even for the yi, which are now of 0{1). In particular, remark that y2 
still converge towards zero, and so does its phase. 

It is thus the larger tuq which should be responsible for the different behaviour shown in 
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Figure 4: RGE evolution of the MFV parameters for the SPS-la benchmark point, but 
with tan/3 = 50. As in Fig. [T], the initial conditions at the GUT scale are varied, but now 
allowing larger range for the ?/», which can be of 0(1)- Still, the quasi fixed-point behaviour 
is obviously largely independent of tan/3. The corresponding plot for 62 is similarly close to 
the tan/3 = 10 one. Those for w\^2 are not shown because these parameters stay very small, 
even with tan/? = 50. 

Fig. El This is confirmed with the help of the analytical solutions, and Table [21 Consider for 
example xi(tew), which now depends ten times more strongly than in the SPS-la case on its 
initial condition xi(to)- If we ignore for a moment that the beta function of x\ also receives 
a contribution from y^ and j/5, which are potentially large, and evaluate the coefficient in 
front of 8\ in the seventh row of Table [21 with the help of Eq. f[331) we get: 



ai(tew) 



h 



'^l(^o) Gnijtc 



6 



(5i0.21 



5.21 
~6~ 



0.18 



(64) 



in perfect agreement with the number in the table. The reason for the enhanced sensitivity 
to the initial conditions is due to the less strong growth of fii in going from the GUT to 
the electroweak scale - the extra dependence on the initial conditions of Xi coming from the 
terms proportional to 1/2 and y"^ is negligible, as we have explicitly checked (the reason can 
be understood rather easily: the i/j's depend on xi through the beta function and vice-versa 
- so the extra dependence of xi(tew) on its initial conditions at the GUT scale comes in 
through the beta function of the beta function and is therefore suppressed by two powers 
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of l/N). In fact, even the value of Xi(tew) for initial conditions of the mSUGRA type can 
be understood rather well on the basis of the analytical formulae obtained in the case of 
moderate tan/?: the relation 2x1 (t) = X2(t) holds to better than one percent and Eq. (I3H1) 
yields numerically: 

xi(tcw) = -(2.10 + 13.76 + 0.42) ■ 10^ = -16.28 • 10^ GeV^ , (65) 

whereas the value obtained solving the exact RGB's numerically is —16.11 ■ lO'' GeV^ - the 
approximate analytical solutions obtained in the moderate tan (3 case work here surprisingly 
well, to one percent accuracy. 

6 Conclusions 

In this paper, we have revisited the formulation of minimal flavour violation (MFV) within 
Minimal Supersymmetric extensions of the Standard Model (MSSM), and linked it to a 
counting rule which keeps track of the hierarchies in the Yukawa couplings and in the CKM 
matrix in a coherent way. This allowed us to move continuously and in a controlled manner, 
between the moderate and the large tan/3 case, keeping control over the expected order of 
magnitude of the different terms. We have argued that to implement these counting rules in 
an efficient way, it is convenient to express the soft SUSY breaking terms of the MSSM in a 
different basis than in the conventional MFV. In order to study the renormalization group 
equations of the MFV parameters, we have projected the beta functions of the soft SUSY 
breaking terms on this basis and checked that the beta functions obey the same counting 
rules as the coefficients themselves. We have then studied the behaviour of the running of 
the MFV parameters numerically, and first checked that this reproduces the full running 
with MFV boundary conditions at the level of accuracy at which we expected them to work. 

In the moderate tan/5 case, we were able to provide approximate analytical solutions to 
the RGB's of the MFV parameters, and we discussed, in the case of the SPS-la benchmark 
point, the behaviour of the MFV parameters under the running. We confirmed the finding of 
Ref. ^ about a quasi fixed-point behaviour of these, and could explain it in detail with the 
help of our analytical solutions. The crucial observation is that the fiavour blind part of the 
soft SUSY-breaking terms, the Oj parameters, run fast, like the gluino masses, whereas all 
the terms with a nontrivial fiavour structure (the Xj's, in our basis) grow less rapidly - their 
beta functions do not contain the gaugino masses, but at best the a^'s. This hierarchy in the 
beta functions explains why the ratios of the fiavour- violating parameters to the fiavour-blind 
ones tend to a small, finite value at the low scale. The only possibility we have identified to 
avoid this, and to make a larger range of values at the low scale at all possible, is to increase 
the value of the Oj's at the GUT scale - this, however, makes the squark masses heavier, 
which in turn suppresses fiavour violations at low energies. 

Our analysis of the SPS-4 benchmark point, which has tan/3 = 50, shows that the picture 
does not change substantially for large tan/?. More parameters (the |/j's in our notation) 
may become of order one, but their RGB's are not particularly different from those of the 
Xj's (they also contain at best the a^'s in their beta functions, and do not grow as fast), nor 
do they influence the other parameters of order one in a substantial way. 
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The picture which emerges from this analysis is that, if MFV is vahd and has its origin 
at the GUT scale (or at any other high scale much higher than the SUSY breaking one), 
then it is a much more constraining framework than if one assumes it to be valid just at the 
electroweak scale. 

This is particularly true for the CP-violating phases. Throughout our study, we have 
allowed MFV coefficients to have imaginary parts, whenever allowed by the hermiticity of 
squark soft SUSY breaking mass terms. The first important observation is that MFV does 
allow for new CP-violating phases not only in the trilinear terms, but also in the LL sector. 
This latter phase could have significant impacts on the MFV predictions for 6 ^ s and b ^ d 
observables (but not for s — > c? ones, where it is absent). However, if MFV has its origin at 
the GUT scale, and tan/3 is not too large, we have shown that all the CP-violating phases 
are strongly suppressed at the low-scale. This behaviour is essentially unaltered when tan (3 
is large, though in that case, the initial conditions at the GUT scale have to be modified. 
In particular, when squark soft SUSY breaking mass terms are as large as in the SPS-4 
benchmark point, the suppression is much less pronounced. 

A final comment about leptons, which in this analysis have only played a marginal role. 
One can of course apply similar ideas in that sector also, as has been shown by Ref. [22] . 
Moreover, if one considers a grand unified theory, one has a more constrained framework, 
and has relations between the MFV parameters in the lepton and in the quark sector [23j 
(indeed, independently from any MFV hypothesis, the relations between the two sectors 
have interesting phenomenological implications, [H]). While these connections are very 
interesting and worth investigating, they mostly concern the boundary conditions at the 
GUT scale - the RGE below the GUT scale are the standard MSSM ones. As a first analysis 
of the running of the MFV parameters, we therefore found it more convenient to concentrate 
ourselves only on the quark sector. 
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A Solving the RGE's and boundary conditions 

For our purpose, only the one-loop RGE's are needed. At that order, the gauge, Yukawa, 
soft-breaking and Higgs sectors essentially decouple, and solving the RGE's proceeds in steps. 
First, the RGE's for the gauge couplings are solved, 

^^ = ^A«- , A = (33/5, 1, -3) . (66) 

Using, the M^S values for simplicity, i.e. a-^ {Mz) = 127.904(17), sin^ Ow {Mz) = 0.23122(15), 
and as (Mz) = 0.1176(20), the unification scale and the corresponding value of the gauge 
coupling are found as 

Mg = 2.0 X 10^^ GeV, a^^Mc) = 24.3 , (67) 

if the MSSM running starts at the Mz scale, which we also assume for simplicity. The 
gaugino masses are required to unify at the same scale Mq as the couplings. 

Ml (Mg) = M2 (Mg) = M3 (Mg) = m,/^ ■ (68) 

Their RGE's are also solved to one-loop. 

Second, the Yukawa couplings at the electroweak scale are set from the known fermion 
masses at that scale (again, we neglect the difference between DR and MS values for sim- 
plicity) 

TTiu {Mz) = 1.27 MeV , m, (M^) = 0.619 GeV , m^ (M^) = 171.7 GeV , 
ma (Mz) = 2.9 MeV , m, (Mz) = 0.055 GeV , rrib (Mz) = 2.89 GeV , (69) 

me (Mz) = 0.4866 MeV , m^ (Mz) = 0.1027 GeV , rrir (Mz) = 1.7462 GeV . 

as well as by setting tan/9 = Vu/vd and t>^ + f J ~ (174 GeV)^. Their one-loop RGE's can 
then immediately be solved numerically. Running up according to either the full one-loop 
MSSM beta functions, or with the approximate MFV RGE's, we find for tan/3 = 10: 

\ / 10^^ 10^^ 10^^ 

I Y^^^ (Mg) I = I 0.00032 0.0014 6.0 ■ 10"^ , AY„ < 10-« 10"^ 10"^ 

0.00459 0.0236 0.577 / \ 10"^ 10"^ 10"^ 

\ / 10^^ 10-9 10-^ 

|Y^^^(Mg)|=| 0.000911 , AYd< 10"^ 10"^ 10"^ 

5.03 • 10-5 0.000259 0.0545 / \ 10"^ 10"^ 10"^ 

Yf^^ (Mg) = diag(0, 0.0040, 0.0675) , AY^ < diag (lO-^ 10-^ 10"^) , (70) 

where AY = |Y^^^ (Mg) \ — |Y^"^^ (Mg) |, and for tan/9 = 50 (neglecting non-holomorphic 
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corrections): 

/ ^ \ ( ^^~^ ^^'^ ^^'^ 

I YMFv ^^^-^ I ^ 0.00034 0.00146 6.2 ■ 10^^ , AY„ < lO^^ lO^^ lO^^ 

\ 0.0050 0.0259 0.683 / \ 10"^ 10"^ 10^^ 

/ \ / 10-^ 10-s 10-^ 

|Y^^^ (Mg) I = 0.00613 , AYd< 10"^ 10"^ 10"^ 

\ 0.00047 0.00242 0.472 / \ 10"^ 10"^ lO^^ 

Y^^^ (Mg) = diag (0, 0.0266, 0.542), AY^ < diag (lO~^ 10"^ 10"^) . (71) 

This shows that our counting is vahd at all scales. It also shows in practice how it is enforced: 
entries of absolute size at or below the O(10~^) are neglected (for consistency, since rriu, rrid 
and me are set to zero), while the non zero entries are reproduced up to (9(10~^) corrections. 
This structure is shared by the squark and slepton soft-breaking terms, following Eq. fl?I]) . 
The evolution of the parameters used to describe the running of the full Yukawa couplings is 
shown in Fig. [5l The scalings of these parameters, in powers of A, can be immediately read 
off these plots. We stress, however, that higher order terms in the RGE, as well as in the 
matching at the electroweak scale, are relevant for the precise values of the Yukawa couplings 
at the GUT scale. In particular, at the one-loop level, the Yukawa couplings fail to unify at 
the GUT scale, as shown in Fig. O This is of no concern for the present work. Rather, our 
purpose was to check whether the beta functions expanded according to our counting rules 
do indeed reproduce very precisely the full one loop RGE. As shown in Eqs. (ITOllTB . this is 
indeed the case, and there is a priori nothing preventing one to extend the present analysis 
to include higher order effects. 

In particular, to show the impact of the non-holomorphic corrections, we consider the 
simulated scenario corresponding to 

A = tan/3(£iYtY„ + £2YjY,) , (72) 

with tan/? = 50 and Si = ±0.002, 62 = ±0.003, which is in the ballpark of the values given 
e.g. in Ref. [21], but still small enough that rji ^ 1 in Eq. 041 p . Such a correction can be 
accounted for by shifting the initial conditions for ys (Mz), yb{Mz) and Cb{Mz). Solving 
Eq. (l38l) iteratively, we find 



£^ = +0.002, 62 = ±0.003 => Vs {Mz) = 0.0138, yb {Mz) = 0.832, c^ (M^) = 0.110 , (73) 
£^ = -0.002, €2 = -0.003 => ys (Mz) = 0.0138, yb (Mz) = 0.683, Cb (Mz) = 0.0595 . (74) 

Taking these initial conditions, and running the Yukawa couplings according to Eq. (I35p . 
or using the full MSSM running, is numerically equivalent to within one part in 10~^, as 
in Eqs. ( ffOl ITTl) . The behavior of the parameters is shown in Fig. O As one can see, the 
parameter Cb(/i), though now slightly larger, is still sufficiently small to be counted as of 
order 0{ybX^). 

The third step is to solve the RGE's for m^^, wi|^^, hiq, m^, m|,, A'^, A^, iii|, m|;, 
and A^, setting initial conditions at the high scale (details of which are given in the text). 
Typically, we start with 

mL (Mg) = mL {Mg) = ml , 
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Figure 5: Up: RGE evolution of the Yukawa couplings, in terms of their skeleton structures 
given in Eg. (1331) . for tan/? = 10 and 50. The running obtained by solving the RGE of Eg. (1351) 
or the full running projected back on Eg. (1331) are indistinguishable. Down: the running of 
the parameters with large tan/3 initial conditions, as specified in Eg. (ff2l) . In the right plot, 
starting with c^ {Mz) < 0, the RGE drives it positive, hence the zero at about 10^^ GeV. 



with ttiq k. Aq also setting the scale of the sguark and slepton soft-breaking terms. 

The final step is to enforce the matching at the electroweak scale for the Higgs sector. 
The Higgs vacuum expectation values Vu and Vd are fixed from tan/3 and Mz- Enforcing the 
correct Higgs potential, and knowing rri^Hu i^z) and rnj^^ (Mz), gives |/i (Mz) \ and b (Mz), 
which can then be run up to the Mq scale. Lowest order approximations are notoriously 
inadeguate at this step, but this is of no concern for us since this matching and running is 
fully decoupled from the rest, in particular from the sguark and slepton sector on which we 
concentrate. 



B Relation between the x^ and the bi parameters 

The standard MFV basis of Ref. pj, and the one with the Xi matrices adopted here are 
related to each other by a linear transformation. The MFV representation of the soft SUSY 
breaking terms given in Eg. ( 12T1) . corresponds to the following in the standard basis (assuming 
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real MFV coefficients for simplicity): 

ai + 6iYtY„ + 62YJY, + ci(yJy,Y1Y„ + YtY^Y^Y^ 
[a2 + 63Y„Yt] , 



m 



Q 



rrin 



m, 



m, 







m 



D 



m, 



A^ 





AoYu 
AoYd 



as 



+ Y, (fee + ferYlY^) Y 

a4 + fegYlY^ + 610 Y^Y^ 

a, + buYlYu + 612 yJy, + ceYjY.YlY^ 



The connection between both bases is explicitly given here: 
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y'^ ci 2cbytybyt 
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\ ctyt Cbyb 2 (qcj, + ybyt) Mt J 
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y^ _ 
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h2 



02 




yhl _ _ 
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VtMb 



-AX^ysym 

/ 1 \ 
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10 

V Q yl J 

Vr yl 

Vii ^ J \ 015 




VbVl 

CbVb 


AX^ysCbyb 

/ Og \ 

bi3 



ctVtyl 

ytytVb 




-AX^ysCbybyl J 



/ «5 \ 

fell 

hi2 



(75) 



(76) 



From this, and the scaling of the Yukawa couplings, one gets the scaling of the MFV co- 
efficients, assuming all a's and 6's are a priori of 0{1). The extension to complex MFV 
coefficients is immediate. 

Note that the trilinear couplings have a non-trivial flavour structure even in the limit 
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bi = 0, and universality is reproduced by setting 




a^Ao 



a^Ao 



agAc 



VtiMc 
CtiMc 
VciMc 
I Vb {Mg 
Cb (Mg) 
Vs [Mg) 




V 




Vr {Mg) 



I 



04^0 



fls^o 
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0.0011 
0.0014 
/ 0.0481 \ 
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V 



04^0 



as^o 



a%A^ 



tan/3 = 50 

0.6393 
0.0052 
0.0015 
/ 0.4163 
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0.5432 
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V 



(77) 



C Higher order terms in the beta functions 

Following the counting rules given in Eqs. (TT71 [22| [3T1 [37|) . the beta functions are expanded in 
powers of A. Actually, this expansion involves only even powers of A, so the first corrections 
arise at (9(A^) when tan/3 ~ A~^. For Yukawa couplings, slepton soft-breaking terms, as 
well as for the Higgs parameters, the beta functions given in the text are already precise to 
(9(A^) or higher, so only those for squark soft-breaking terms are missing. Denoting b(5c the 
C(A^) corrections, they are easily found to be 

5/3al = 5/3a2 = <5/?a3 = , 

5/3^1 = 2 (cfei/6 + CtVt) {V\ + Rey2) , 
^i^yi = 2 (cbyb + ctVt) iyi + Re 1/2) , 
5(3y^ = {ctyt + Cbyb) (2ai +xi+yi + 2y2) + 2ctyt (m^^ + 02 + X2) + 2cbyb {m]j^ + 03 + I/3) , 

%2 = ^CtVt {mjj^ + di + a2 + xi + X2 + yi + 2 Re 1/2) 

S(3y^ = Scbyb {m]j^ + ai + 03 + xi + yi + 2 Re 1/2 + Z/3) , 

Spy,, = -2A\^Cbys {2m]j^ + 2ai + 2a3 + 2xi + 2y2 + 2/3) - 2AA^ ((xi + ^2) ybVs + w^y^) , 

5(5a^ = tti {25ctyt + Cbyb) + y^ {Hctyt + Cbyb) + 2y5 {ctyb + Cbyt) , 

5j3y^ = 04 (llctyt + Cbyb) + 1/4 (25Qyi + Cbyb) + 2a^ {ctyb + Cbyt) + CtK'^ , 

5f5w^ = 6ct ((04 + 2/4) yc + W2yt) , 

%5 = «5 {2hcbyb + Ctyt) + y^ (Hcfeyf, + ctyt) + 2yi {ctyb + Cbyt) , 

5(3y^ = ar, {llcbyb + ctVt) + 2/5 {25cbyb + ctyt) + 2a4 {ctyb + Cbyt) + Cb {2asyr + K'J , 

5^^3 = Qcb ((as + 2/5) ys + w^^yb) , 

5(3w^ = 2w4 {5cbyb + ctyt) - AX^ {{5cbyb + ctyt) W3 + Acbys (as + Vb)) 



2p + 2ir]) yt {w^yt + 2 (a4 + 2/4) ys) ■ 

(78) 
Including these corrections, the beta functions are all correct up to tiny 0{\^) corrections. 
Numerically, they reproduce the full MSSM one-loop running, projected back on the Xj basis. 



34 



to better than 1%, and even to better than 0.1% for mSUGRA-type of initial conditions at 
the GUT scale. For the numerical analysis discussed in the text, we have always used the beta 
functions including the corrections given here - we stress that even if formally suppressed, 
the terms containing K'^ or m^j may be numerically important, whereas most of the others 
give very small contributions. 

D Fixed points 

The numerical analysis has shown that all the MFV parameters tend to a certain "fixed 
point" at the low scale - rather strongly for the SPS-la point and less so for the SPS-4 one. 
Here, we discuss this feature in more detail, and try to gain some analytical understanding of 
it on the basis of the simple beta functions derived in Sect. [SI A true fixed point occurs if the 
beta function of a parameter has a zero which depends only on the value of that parameter 
- depending on the sign of the derivative of the beta function with respect to the parameter 
at the position of the zero, the fixed point is an ultraviolet (positive) or an infrared one 
(negative). In our case, we observe the behaviour of a quasi infrared fixed point. 

Such behaviour is observed for the ratios Xi/oi and X2/a2, for example, but also for all 
other similar ratios. The beta functions for these ratios read: 



Hxi/ai — ~ 



Pxi ~ PcLi 



(79) 



and we have to find out whether these beta functions have zeros, and on what parameters 
these zeros depend. The equations [3x^/ai = are unfortunately nonlinear, and all coupled 
to each other, such that an analytical solution is difficult to obtain in general, and in any 
case, not very illuminating because too cumbersome. The terms which make the equations 
nonlinear, however, all originate from the dependence of the Sj's parameters on the Xj's, etc., 
which, as we have seen, is negligible. The behaviour of the Sj's is mostly driven by the gluino 
masses, and somewhat also by the terms proportional to yt. In this approximation, it is then 
easy to solve the equations and obtain simple analytical expressions for the solutions. In 
addition, if one considers only moderate tan/3 and neglects higher orders in A, the equations 
simplify further, and look as follows: 

(80) 

Analogous results can be given for all the other parameters. These zeros of the beta functions 
are not true fixed points, because they depend on parameters which do run - they are rather 
some sort of "running fixed points" . This means that even if a parameter reaches exactly its 
"fixed point" a certain scale, it will not stay there because the zero of the beta function will 
move with the scale. On the other hand, even if moving with the scale, they do represent a 
line of attraction for the parameters. 

For the ratios xi/fii and X2/a2, the running fixed points fiatten out at the low scale, 
such that in that region they become almost true fixed points. Indeed, also the numerical 
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Xi 


"i^Hu +0.1 + 0,2 + o^/y^ 


X2 _ 


-2^' 


ai 


~a, + 2^2 + ^glMllyl 


02 


ai 



evaluation shows that they almost coincide with the values these parameters tend to (for the 
SPS-la point): 

^ = -0.18 , ^ = -0.18 , 

-± = -0.39 , :r^ = -0.38 . (81) 

«2 \mz «2 Im^ 

These formulae do not allow one to understand how strong these quasi fixed points are at 
the low scale - the proper explanation of the behaviour of the RGE has been provided in 
Sect. O 
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